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Abstract. A difference equation analogue of the Knizhnik - 
Zamolodchikov equation is exhibited by developing a theory 
of the generating function H(z) of Hurwitz polyzeta functions 
to parallel that of the polylogarithms. By emulating the role of 
the KZ equation as a connection on a suitable bundle, a differ- 
ence equation version of the notion of connection is developed for 
which H(z) is a flat section. Solving a family of difference equa- 
tions satisfied by the Hurwitz polyzetas leads to the normalized 
multiple Bernoulli polynomials (NMBPs) as the counterpart to 
the Hurwitz polyzeta functions, at tuples of negative integers. A 
generating function for these polynomials satisfies a similar differ- 
ence equation to that of H(z), but in contrast to the fact that said 
polynomials have rational coefficients, the algebraic independence 
of the Hurwitz polyzeta functions is proven. The values of the 
NMBPs at z = 1 provide a regularization of the multiple zeta val- 
ues at tuples of negative integers, which is shown to agree with the 
regularization given in |AET01j . Various elementary properties of 
these values are proven. 
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The equation VG = satisfied by the flat sections of the universal 
prounipotent bundle XL with connection V on P 1 \ {0, 1, oo} is known 
as the formal Knizhnik-Zamolodchikov (KZ) equation. A distin- 
guished solution is given by the polylogarithm generating series Li (2), 
which is the Chen series 
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where the integrals are iterated, the sum is over all words W in non- 
commuting variables X and X\, and if u and lo\ denote — and 
respectively, for any word W — X h . . . X ir then LOW '. — . . . USi r . The 
KZ equation arises in the study of the shuffle algebra of polyzeta values 
- see |Car01] for the details. 

The goal of the present note is threefold: to introduce an analogue of 
the KZ equation which emerges from purely algebraic considerations 
and is related to the stuffle algebra of polyzeta values, to develop the 
notion of difference connection modelled on this equation to mimic the 
usual idea of (differential) connection on a manifold as in the case of 
the KZ equation, and finally to discuss normalized multiple Bernoulli 
polynomials in this context. Here are some details: 

The role of the polylogarithm functions in this setting is played by the 
Hurwitz polyzeta functions 



where the sums converge for (sj + . . . + s r ) > r—j + 1 and determine 
meromorphic functions of z with poles at non-positive integers. The 
Hurwitz polyzeta functions may be seen to give rise to a homomor- 
phism from a stuffle algebra of polynomials over C in infinitely many 
non-commuting variables {yj}jZi, to a stuffle algebra of meromorphic 
functions on C - i.e. 

where C(Vji ■ ■ ■ Vj r \ z ) * s some (sum of) Hurwitz polyzeta function(s). 

Introducing dual variables {Yj}'^ =1 to the yj, the algebra of non-commuting 
power series in the Yj with some ring A of coefficients may be endowed 
with the structure of Hopf algebra. When A comprises a class of mero- 
morphic functions of which the Hurwitz polyzetas form a subclass, a 
distinguished element of this algebra A is the generating series 



where the sum is taken over all words in the yj. From formal consid- 
erations (cf. |Rac00j ). it is evident that H(z) is group-like - i.e. with 
respect to the comultiplication A in A, 




h(z)-.= c(v*...v*i*)iSi 



Vil-Vjr 



AH(z) = H{z)®H{z). 
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Moreover, making use of a family of elementary difference equations 
satisfied by the Hurwitz polyzeta functions, we prove the 

Theorem A. If D_ denotes the difference operator D_F(z) = F{z) — 
F(z - 1), then 

OO y 

(!) ^-^H-Et^zW^)- 

k=l 1 ' 

We go on to define an algebraic analogue of the notion of connection, 
for which <^ determines some kind of universal algebraic connection 
V alg with respect to which H(z) forms a flat section. Emulating the 
unipotence property of connections in this context, and defining an %R Q - 
bundle over the complex sphere with finitely many points removed to 
be any bundle with the space of meromorphic functions on the sphere 
with poles exactly at the non-positive integers as subbundle, the versal 
property is the content of the 

Theorem B. Given any DJlo-bundle £ with unipotent difference con- 
nection V al9 on X along with a point b 6 X and any v 6 there 
exists a mapping 

^:(U,V A )^(£,V^) 
which is compatible with the difference connections and has (i/j v )i,(1) = 
v. 



In the final section of the paper, a generating function for normal- 
ized multiple Bernoulli polynomials is defined, and shown to satisfy 
a similar difference equation to ([I]). Precisely, the above-mentioned 
difference equations may be considered for negative arguments. Us- 
ing known solutions for the equations of first level given by the usual 
Bernoulli polynomials, these difference equations may be solved ex- 
plicitly. The polynomials ((— ni, . . . , — n r \z) which arise are counter- 
parts at negative integer parameters (— ni, . . . , — n r ) to the Hurwitz 
polyzeta functions defined above. Denoting by {Yj}j< a countable set 
of non-commuting formal variables, one then sets 

H B (z) := ^ C(-Wi> • • • > -n r \z)Y_ ni . . . Y_ rir 

where the sum is taken over all tuples (— n 1; . . . , —n r ) of non-positive 
integers. As before it is possible to prove 

Theorem C. 

oo 

D_H B (z) =-J2 Y- k (z - l) k H B (z). 

k=0 
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It is interesting that while these normalized Bernoulli polynomials 
have rational coefficients, the Hurwitz polyzeta functions are tran- 
scendental. Indeed, we prove: 

Theorem D. The Hurwitz polyzeta functions are algebraically inde- 
pendent over C - i.e. for any N > 1, if P(T%, . . . ,T N ) G C[T 1; . . . ,T N ] 
is a polynomial which vanishes (uniformly in z) at some N -tuple of 
any Hurwitz polyzeta functions of the variable z, then P is identi- 
cally zero. 

Evaluating the Hurwitz- Bernoulli polyzeta functions at z — 1 gives 
the polyzeta values (also known as multiple zeta values in the litera- 
ture) in the case of the Hurwitz polyzetas, and a set of possible reg- 
ularizations of the polyzeta values at non-positive integers in the case 
of the normalized multiple Bernoulli polynomials. It is known that 
many different regularizations are possible - cf. |GZ08j - and in the 
case of our difference equation method, it is clear that the values are 
not uniquely determined since the constant term of any polynomial 
vanishes under the action of DA However, it is possible to show that 
these regularizations coincide with those given by a different method 
in |AET01] . More precisely, by determining a recursive formula for the 
normalized Bernoulli polynomials, one can prove the 

Theorem E. 

£(— k\, . . . , — k r \l) = lim ••• lim C( s i> ■ ■ ■ > s r)> 

Si— > — ki s r — >— k r 

where C(si, . . . , s r ) denotes the meromorphic continuation of the series 



Once this is proven, we establish various properties of the regular- 
ized polyzeta values using the explicit normalized multiple Bernoulli 
polynomials, for example 

Proposition A. If n and k are integers for which n > 1 and < k < 
2n + 1, then 




C(0)C(-2n - 1) = C(-2n + k — l, -k\l) 



and 



Proposition B. If n > 1 is any integer, then 

C (_ n) _ 7l _l|l) = C (_ n _l > _ n |l). 
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1. The Hurwitz polyzeta functions and the formal KZ 

DIFFERENCE EQUATION 



With notation as in the introduction, on the space A << yi, y 2 , . . . >> 
of formal power series in the non-commuting variables {yj} c *L l as in 
[CarOl] define the stuffle product by induction via 

y k w * y k ,w' = y k {w* y k >w') + y k >(y k w * w') + y k+k >{w * w'), 

where 1 denotes the empty word, linearly extended to arbitrary for- 
mal power series. Those words in the yj which do not end in yi may 
be mapped to Hurwitz polyzeta functions via the homomorphism 
described above. For words which end in yi, some regularization is 
required. This is achieved by setting 

sv 1 ; r(z) 

where T(z) denotes the usual interpolation of the factorial function; 
then building up the remaining functions for which some regulariza- 
tion is required, using the stuffle product itself. The choice of this 
regularization is motivated by the classical expression: 

1 T'(z) 

^ = — -tU + 0{ '- 11 

i.e. 

(2) .im(W;,)-i)=-^ 

This is an apt choice for our purposes since Clll- 2 ) as defined here 
satisfies a functional equation of the form of similar equations satisfied 
by the other functions in the family of Hurwitz polyzeta functions: 
Directly from the definitions one may deduce: 

(3) C(h, K\z + 1) - C(h, ...,k r \z) = — j-C(A&, ...,k r \z+l), 



where if r = 1, 



(4) C ( fc | 2 + i)_ C (A;| 2 ) = _^ forfc>l. 

z 

An elementary calculation using the functional equation T(z + 1' 
zl(z) shows that also 

c(i\z+i)-ai\z) = ~. 
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Remark 1.1. The regularization procedure discussed above admits of 
an interesting formulation involving concrete limits in the style of d2J): 
To see this, we introduce non-commuting variables y t for real t > 1 
and define a stuffle product on the non-commuting polynomial algebra 
A < {yt}t>i > as before: 

Vh * yt 2 = VhVh + yt 2 yh + yt 1+ t 2 

for any t\\ t 2 > 1. We denote this algebra by f) s t,R. Also, write 

lim 7/ fc+£ = y fc 

for any k G N\{0}. Such limits clearly commute with the stuffle prod- 
uct. 

Now consider the Hurwitz polyzeta functions C(si, . . . , s r \z) where 
the Sj G R have Sj > 1 for j = l,2,...,r — 1 and s r > 1. Such 
functions satisfy the stuffle relations and form a Q-algebra which is 
isomorphic to a subalgebra, say f)!? t R , of P) s t,R ? when we take A = Q. To 
extend the correspondence to the entire f) s t,R, as before we employ the 
regularization for C(l| z ) given above. Then (j2J) is the analogue for the 
regularized Hurwitz polyzeta functions of 

lim y 1+e = yt 

and evidently commutes with stuffle product among these polyzeta 
functions. In this way, any regularized Hurwitz polyzeta function 
(at a tuple of integers) may be expressed as a limit which equals the 
stuffle product expression that can be built up as before to give the 
regularized value. The point is that regularization assigns meaning to 
a certain symbol, (which amounts to the assignment of a function of z in 
the instances of the Hurwitz polyzeta algebras), and here, via certain 
well-defined limits an alternative description of the regularization is 
supplied. An example will clarify these ideas: 

Example 1.2. Here the limit regularization of £(k, l\z) for any integer 
k > 1 will be elucidated: Corresponding to the stuffle product 



Vk * Vl+e = Vkyi+e + Vl+eVk + Vk+l+e 



AN ALGEBRAIC KZ EQUATION 



we have the equality 
(5) 

E 



+ E 



+E 



(z + n 1 ) k (z + n 2 ) 1+e ^ (z + n 1 ) 1 + £ (z + n 2 ) k ^(2 + 

0<ni<n 2 V ' V ' 0<ni<n 2 V ' V 1 n=0 K 

for e > 0. 
Then 



E — I IE 1 ) - (E — I (- 



[z + n) k I \ ^ (z + n) 1+£ I \^ (z + n) k J \e 

v n=0 v ' / \n=0 K ' / \n=0 V ' ' ' 



1 



and in the limit as e — >• 0, this expression approaches 

r(zY 



y— 

^ (z + n) k 



,n=0 



Viz) 



ak\z)ai\z) 



Solving (JSJ) for the sum we are attempting to regularize, namely 

1 



E 

0<ni<n2 



(z + n 1 ) k (z + n 2 



|1+E 



and subtracting 



oo 1 1 



n=0 



from both sides of the result before finally taking the limit, we find 
that 



C(k\z] 



lim C(k,l + e\z) - 



C(fck)C(ik)-C(i,fck)-C(* + i|«), 



which is the regularized value of l\z) 
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As a sampling of similar computations we mention that for any integers 
k\ , . . . , k r where k r > 1 , 

lim (C(h, ...,k r ,l + e\z)- C(ku...,K\z) \ = _ k r \z)((l\z) 

r r 

— ^ ] C(^-i) ■ ■ ■ ) %— i) I; kj, . . . , k r \z) — ^ ] C(ki, ■ ■ ■ , kj-i, kj + 1, fcj+i, . . . , k r \z). 
i=i i=i 

Also, in regularizing £(2, 1, l\z) we obtain 



lim lim C(2, 1 + e, 1 + <$U) + C(2, l + S,l + e\z) 



5\ £ J 5 e 5 e 

= C(lk) 2 C(2k) + C(l|*)(-2C(3|z) - 2C(1, 2\z)) + C(4|z) + 2C(1, 3|z) + 2C(1, 1, 2|«). 

1.1. The algebraic KZ equation. Introducing variables {1^ }^ =1 dual 
to the yj, we may form the graded Hopf algebra dual to the algebra 
A < jjj >j*Li of polynomials in the non-commuting yj. In this case, the 
comultiplication is given by 

fc-i 

A ST (n) = l®Y k + J2 Y i® Y k-j +Y k ®l. 

3=0 

This is designed in such a way that in each term on the right hand side, 
the stuffle product of the two constituents of the tensor product has 
the left hand side as a term. As a consequence of this it is not hard to 
show by a formal argument that the generating function of Hurwitz 
polyzeta functions 

H ( z )'-= C(yk 1 ---yk r \z)Y kl ...Y kr , 

(where the coefficient of the empty word is set to be 1), is group-like - 
i.e. satisfies 

A ST H(z) = H(z)®H(z). 

The analogue for H(z) of the differential equation satisfied by Li(z) 
is, as disclosed above, of an algebraic nature: The equations (EJ) can 
be brought together in a universal functional equation for Hurwitz 
polyzeta functions, just as the formal KZ equation is a compilation of 
differential equations satisfied by the various polylogarithm functions. 
To be precise, we have the 
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Theorem 1.3. If D- denotes the difference operator D -F (z) = F(z) — 
F(z - 1), then 

oo 

(6) D-H{z) = -^j^H{z). 

k=i 1 ' 

Proof: From ([3]), it is clear that for fixed k > 1, we have 

C{y k w\z)Y k w = ^ w \ z ~ l ) Y k w ~ T^rW E ch*) 
= J2((yMz-mw--^H(z) 

where W is the word in the variables Yj which is dual to w. 
But then certainly 

oo 

k=i ^ ' 

and the statement of the theorem is immediate. □ 



2. Difference connections 

The Knizhnik-Zamolodchikov equation characterizes flat sections 
of the universal unipotent bundle with connection on P 1 \ {0, 1, oo}. 
Although much less rich, an algebraic analogue can be ideated around 
(JOJ) - To this end, use DJl k to denote the family of meromorphic functions 
of a single complex variable which have poles at all integers less than 
or equal to some integer k. Then T(z) G Tt as is ((si, . . . , s r ; z) 
(for suitable tuples of complex Sj) when viewed as a function of z. 
Notice that if F(z) G VJt k , then F(z — 1) G SDtk+i, so the operator 
D~ : F(z) i — y F(z - 1) sends ffl k onto Tl k+1 . Now let Wl^ denote the 
union of the DJl k for all integer k, a space which is invariant under D~. 

We have D~ = I — D- and ([6]) is the same as 

"-'-E^) *«=<>■ 

The operator D~ — I — Y^k=i (z~i) k s ' lou ld play the role of an algebraic 
version of a connection. 

An OJtoo-bundle £ on X := Pj.\{a 1 , . . . ,a N } where the aj G Z U {oo} 
is a bundle on X having WX k as a subbundle for each integer k, where 
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by an abuse of notation Wlk =■ Wlk(X) comprises those functions on 
X which when viewed as functions on C using the usual coordinates 
are meromorphic with poles at all integers less than or equal to k. (We 
also require a pole at oo if the latter is a point of X.) An 9Jt fc -bundle 
on X is a bundle merely required to have DJlk itself as subbundle. 

Notice that QJtoo =: VJt^X) may itself be regarded as an DJtoo-bundle 
on such an X. This bundle is the analogue of the space of differential 
1-forms on X. 

With this notation, we make the 

Definition 2.1. Given an 9Jto-bundle £ on X, a difference connection 
on the bundle is a C-linear mapping 

V al9 : £ -> OJtoo ® OTo £ 

for which for any / G $Jloo and any section s of £, 

(7) V al9 (fs) = D-f ■ (-Ds) + f(V al9 s). 



The condition (jTJ) should be regarded as an algebraic analogue of the 
usual Leibnitz rule. 

In what follows, any bundle will be assumed to be an 9Jto-bundle. 
We require the following notion: 

Definition 2.2. A bundle with difference connection (£, V a ' 9 ) is called 
unipotent when there is some r and some sequence of upper triangular 
matrices {Nk}k>i such that 

°° N 

(£,V^)^«^--/-Et7^)- 



We proceed to construct an object which satisfies a versal property 
with respect to the bundles with unipotent difference connection. 

Consider the algebra C < Y > of polynomials in the countable set of 
non-commuting variables Y := {Yj} c *L 1 . Let J denote the augmentation 
ideal J = (Yi, Y 2 , . . .). Then write 

U n := C < Y > /J n+1 . 

This algebra comprises those polynomials in which the words that ap- 
pear have length at most n. Now form U n := f/ n ®9Jto an d U = lim^ U n . 
Observe that we could write U = WIq «Y» since effectively we are 
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considering the formal power series algebra in the non-commuting vari- 
ables Y k with coefficients taken from 9Jt . Also, notice that H(z) G U 
when we take X to be C. 

On U a difference connection V" 4 is determined by defining for each n, 

\ \w \<n 

|iu|<n |w|<n \w\<n k=l 

where [w] denotes the class of the word w, and pr n indicates projection 
to words of length less than or equal to n. 

Theorem 2.4. Given any DJl -bundle £ with unipotent difference con- 
nection V alg on X along with a point b G X and any v G £;,, there 
exists a mapping 

ij v : (U,V A ) ->■ (£,V a/s ) 
which is compatible with the difference connections and has (ip v )b(l) = 
v. 

Proof: Corresponding to a word w = Y h . . .Y ir , set N w := N h . . . N ir . 
Then define 

i) : (U, V A ) -> (£,V ais ) 

by taking 

\ TO /to 

Notice that = 1 • t> . 

Now N k N w = Ny k w Because of this, one sees readily that ip o V A = 
o \7 alg . □ 

As in the topological case, one can define a unipotent fundamental 
group n^ R ' dl ^ (X, b) as the tensor compatible automorphisms of the 
fiber functor with respect to the category of unipotent difference con- 
nections. This group acts on sections of any bundle with difference 
connection, and by definition, the action commutes with the ifi v of the 
Theorem. This gives the "parallel transport" action on the bundle, but 
unlike the iterated integral situation, this action admits no description 
intrinsic to the flat section of the bundle (with respect to the difference 
connection). Moreover, while the parallel transport provides a means 
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of showing that the analogous mapping to ip v in the topological con- 
text is characterized by its action on the element 1 of the fiber above 
b, via the linear independence of the iterated integral (polylogarithm) 
functions, even though the Hurwitz polyzeta functions are linearly in- 
dependent over C, this linear independence does not interact suitably 
with the action of k 1 ' (X, b) to facilitate similar conclusions. 

By construction, H(z) is a flat section of (IX, V A ). 



3. The normalized multiple Bernoulli polynomials 
The classical Hurwitz zeta function defined for IHe s > 1 by 



oo 1 



, .z + n) s 

?1=0 V ' 

may be analytically continued in the complex variable s: For example 
modifying Riem Ann's contour integral approach to his zeta function 
C(s), it is easy to show that 

r(l-s) f eQ--*> dw 



27ri J c e w — 1 w 

where C is the Hankel contour in C (i.e. a loop about based at 
infinity, enclosing the positive real axis). This integral expression con- 
verges for all values of s G C\{1}, and (as was known classically) at 
non-positive integers s = —k one finds 

(8) C( _ M = _:?5±iW 

where Bk(z) is the fcth Bernoulli polynomial. (Again cf. [CarOlj .) 
Henceforth write ((— k,z) =: ((—k\z). 

These normalized Bernoulli polynomials thus belong to the family of 
Hurwitz zeta functions, and if ((s) denotes the Riemann zeta func- 
tion and B k := B k (0) = (-l) k B k (l) is the fcth Bernoulli number, 
then the fact that 



C(-*) = (-i 



k+ 1 

concords with ((s, 1) = C( s )- 

Various generalizations of these polynomials to multiple versions exist, 
going back to work of Barnes in 1899 and much more recently Szenes 
in |Sze98j . and Komori, Matsumoto and Tsumura in |KMT10j . 
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Our approach here is somewhat different however, being determined 

by solving the system of difference equations 

(9) 

V(-k u -k r \z+l)-V (-kx, -k r \z) = -z kl V(-k 2 , -k r \z+l) 

where the kj are non-negative, for all possible r > 1. Of course, these 
are the difference equations (j3J) satisfied by the Hurwitz polyzeta func- 
tions, but at parameters (—hi, . . . , —k r ) at which the polyzeta functions 
do not exist, (cf. [AKT01] ). 

Theorem 3.1. Up to addition of a function of period 1, the solutions 
to are given by the recursive formula 

C( ki, . . . , fc r |z) = - ■ ~C( k\, . . . , k r —2, k r — % k r 1|^) 

k r + 1 

1 ^ r 

-K-2, -k r -i-k r \z)+^2(-k r )g -. ~T\j C(-fei ; • • > -K-2, -k r „i-k r +q\z) 

q=l ^ >' 

writing ((—k\z) = — Bk £+^ as above. 

Proof. As is shown in |Mes59j the (normalized) Bernoulli polynomi- 
als £(— n\z) give the unique solution to 

(10) u(z + 1) - u(z) = -z n 

for any n > 0, up to addition of a periodic function - i.e. the solution 
to (E]) when r = 1 corresponds to £(— k\\z) = — Bk ^+^ ■ 

Since each B n (z) is polynomial, when r = 2 the right side of (Q is 
also polynomial. Consequently, by linearity, (j5J) reduces to a sum of 
equations of the form of (fTUj) . so that some linear combination of usual 
Bernoulli polynomials gives a solution. Then for r = 3, again the 
right hand side of (|9]) is polynomial so may be solved by the same 
method. Proceeding inductively, it is clear that polynomial solutions 
exist for all possible r, which may be explicitly determined: 

Consider firstly the case that r = 2, and recall the well-known formula 



3=0 V 7 

Using ( ITU]) with n = k 2 , © then becomes 
V(-h, -k 2 \z+l)-V(-k 1 , -k 2 \z) = z ^+ 1 -±- [ ( h + 1 ) BjZ 



ki+k2+l—j 



14 



SHELDON JOYNER 



for which a solution is 

1 k2+1 ( k 4- 1 \ 



Now write out the j = and j '• = 1 terms of the sum, set / = j — 1, 
and use the Pochhammer symbol notation to write 

k 2 (k 2 -l)...(k 2 -(l-l)) = (-k 2 ) l (-l) 1 . 

One finds 

C(-h, -k 2 \z) 

= ~C(-h - k 2 \z) - j^-rCi-h ~k 2 - l\z) 
2 k 2 + 1 

+ • (^yyy ■ C(-*i - A* + 

The (— 1)' +1 factor of the last line may be omitted since -B;+i = for 
all even I > 2. 

Continue by induction, supposing that it is known that a solution to 
([9]) has the requisite form for all positive integers r less than or equal 
to some N — 1. Then by this inductive hypothesis, 

C(-h, -k N \z + 1) - C(-fci, • • • , ~k N \z) 
= -z k X(-k 2 ,...,-k N \z+l) 
1 



-z kl 



k N + 1 



C(-k 2 , . . . , -fc r -2, -fcjv-i - k N - l\z + 1) 



1 

- 2^{~k 2 , -k N - 2 , -k N -! - k N \z + 1) 

k N n 

+ ^(-fciv)g ^ ^Tj? C(-^2, • • • , -^7v_2, -fcw-i - ^iv + ?k + 1) 

The difference equation obtained here may be regarded as a sum of 
difference equations of the form of 

v(z + 1) - v(z) = -Az kl w(z + 1) 

where in each case, w(z) is some normalized multiple Bernoulli poly- 
nomial of depth N — 2 and A is some rational number. Each of these 
may be solved and the sum of these solutions gives the desired formula 
in the case that r = N. □ 
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Corollary 3.2. 

({-fa, . . . , -k r \l) = lim ••• lim ({s h . . . , s r ), 

si— ki s r — ¥— k r 

where ({s\, ■ ■ ■ ,s r ) denotes the meromorphic continuation of the series 

1 



n-t . . . ni r 

0<m...<n r 1 r 

This is immediate from the Theorem combined with the expression for 
the limit given in |AET01j . 

3.1. Elementary facts pertaining to values of NMBPs. 
Proposition 3.3. 

C(0)C(-2n - 1) = C(-2n + k — l, -k\l) 
for n > 1 and < k < 2n + 1. 



Proof. 



2n + 2 fc+l^V ^ / 2n - ? + 3 



so that 



3=0 v 7 J 



In the sum that appears on the right side, since j and 2n — j + 3 have 
opposite parity, the only non-zero term is that for which j = 1. (Notice 
that 2n - j + 3 > 2.) 

Hence, 

C(-2rH-*-l, = ^±A^± = _I = C ( )C(-2n-l). 

sv ' Iy 2n + 2 22n + 2 2 V 2n + 2 J sv AV ; 

□ 

Proposition 3.4. For any integer n > 0, 

C (_ nj _ 7l _l|l) = C (_ n _l > _ n |l). 
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Proof. Using 

1 / n + 2 \ 1 / n + 1 \ ( n \ j — 1 



n + 2\ j y n+l\ 3 J \3 J {n + 2 - j)(n + l - j) 
one computes 

C(-n, -n -l\z)- C(-n - 1, -n|*) 



ti\ j-1 D B 2n+ 3-j(z) B n+2 B n+1 (z) 

■JD-j ■ — r + 



n+l 

j J (n + 2-j)(n+l-j)^ 3 2n + 3 - j ' n + 2 n+l 

Now evaluate at z = 1 and observe that the parity of j and 2n + 3 — j 
differs, while 2n + 3 — j > 2 since n > 0. Also, the parity of n + 2 and 
n+l differs. Thus we can use the same argument as in the proof of 
the previous proposition, to conclude that 

£(_ n> - n - 1|1) - C(-n - 1, -n\l) = 0. 

□ 

Proposition 3.5. 

C(-2m - 1, 0, -2n 2 — 1 1 1) = -C(-2ni - 1, -2n 2 - 1|1) 
/or any ni, n 2 > 0. 

Proof. Once the polynomials ((— 2n\ — 1, 0, — 2n 2 — l|z) and £(— 2n x — 
1, — 2n 2 — l|z) have been determined explicitly, evaluation at z = 1 
gives an expression which may be readily simplified by considering the 
parity of the indices j, k in the products Bj ■ B^ which arise, exploiting 
the fact that the only odd j for which Bj is non-zero, is j = 1. One 
obtains 



C(-2m-l,0,-2n 2 -l|l) 



2n 2 +2 



B 2ni +2n 2 +3 



2ni + 2n 2 + 3 2n 2 + 2 ^— ' V m / m 2n l + 2n 2 + 4 - m 

1 z z m=0 ' 

-C(-2ni - l,-2n 2 - 1|1). 



□ 



It is likely that also 

C(-2n + k, 0, -k\z) = ({-2n + k, -k\z) 
for any k, n with < k < 2n. 



AN ALGEBRAIC KZ EQUATION 
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It is convenient to introduce the following notation: If f(z) = Yln=o a n zn ■, 
write f{B)(z) for Ylt=o a n Bn ^[ z) ■ Also, for ((m, . . . , n r \z) with m = 
. . . = n r = n, write ( r (n\z). 

Proposition 3.6. 

Cr(0|l) = (-I)" ' 



r + 1 
and 

Cr(0|0) = (-l) r+1 ^TT- 

The first assertion is proven (by different means) for the limit expression 
of Corollary S3 in jATOlj . 

Proof. First observe that by applying the difference equation satisfied 
by the multiple Bernoulli polynomials recursively, 

r-1 

( r (0\z + 1) - Cr(0|*) = X;(-l)"Cr-»(0|z). 

n=l 

Hence 

r- 1 

Cr(0|z) = X;(-l) n Cr-»(0|B)(«). 

n=l 

But then taking t = r and t — 1 = r in this equation successively and 
adding, 

(ii) CMz) + Ct-im = -Ct-MB){z). 

Now since C(0|<z) = — z + |, 

C(o,o|z) 



^ 2 l 



2 6 
and 

z 3 z 2 z 1 
C(0,0,0|z) = --- T + - + -, 

the assertion is seen to hold for r = 1,2,3. Now when the statement 
holds for both r and r + 1, then by ( ITT]) . 

(-l) r+1 2 
r + 2)(r + l)r' 

which is true in the case of r = 1, 2 by the above. 



(12) -Cr(0|S)(0) 
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We proceed by induction, supposing that ( 1121) is known to hold for all 
r less than or equal to some K. Also suppose that the assertion of the 
theorem is known to hold for such r. 

Now ( ITTj) holds for all t - in particular for t = K + 1. But then 

a+i(o|o) = -Cjf(o|o)-CA-(o|s)(o) 

{-l) K+l {-l) K+l 2 
~ ~ K(K + 1) + (K + 2)(K + 1)K 

(-l)K+2 



(K + 2)(K + iy 

proving the assertion regarding evaluation of £ r (0|2) at z = 0, by in- 
duction. 

For the other statement, recall the difference equation: 

( K+1 (0\z + 1) - Ck+x(0\z) = -<Zk(0\z + !)• 
Evaluating at z = gives 

a+i(o|i) = a+i(o|o)-a(o|i) 

= — as above, and by induction hypothesis 

(K + 2)(K+1) K + l 

K + 2 ' 

□ 

Remark 3.7. One can easily deduce from the above proof that also 

f-l) r+1 2 



-Cr(0|£?)(0) = Cr+l( 



' (r + 2)(r + l)r 
for r > 1. 

Proposition 3.8. If (—n\, . . . , — n r ) is a tuple of non-positive integers, 
and —ni < 0, then 

C(-ni, . . . , -n r \l) = C(-m, • • • , -n r |0). 

Proof. Because -02n+i = unless n = 0, the assertion is known from 
the equality B n (0) = (—l) n B n (l) in the case that r = 1. 



The general statement now follows from the formula in Theorem 13.11 
by an easy induction. □ 
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3.2. The algebraic Knizhnik-Zamolodchikov equation for MBPs. 

Let {?/_ n }^L denote an infinite family of formal non-commuting vari- 
ables. To any monomial in the jjj we associate a normalized multi- 
ple Bernoulli polynomial (NMBP subsequently) in the obvious way, 
namely 

Vh ■ ■•Vir *-> CiVh ■ --VM '■= C(~h, • • • , -i r \z). 

Now take {Y_ n }^L to denote some other family of formal non-commuting 
variables. Then form the generating function of NMBPs: 

H B (z) := C{w\z)W 
weY* 

where the sum is taken over all words w in the yj and W denotes 
the corresponding word in the Yj (so that ■ ■ ■ y ir is associated to 

Theorem 3.9. 

oo 

H B {z + 1) - H B {z) = -J2 Y-kZ k H B {z + 1). 

fc=0 

Proof: Fix k > 0. Then 

£(y-kw\z+l)Y- k W- C(y~ k w\z)Y. k W = -z k Y„ k £{w\z+l)W. 
weY* w&* weY* 

Adding over all k > gives the result of the theorem. □ 

This proof is virtually identical to that of Theorem 11.31 for the corre- 
sponding fact for Hurwitz polyzeta functions. Consequently, the next 
result follows trivially: 

Corollary 3.10. 

oo 
k=0 

is a difference connection on P 1 . 

Now it is very tempting to try to form a generating function encom- 
passing both classes of functions discussed so far. However, this would 
necessarily include polyzeta functions of the form of C(^i> ■ ■ ■ ,k r \z) for 
arbitrary tuples of integers kj, not all of which can be determined by 
difference equation methods. In fact, by a similar recursive procedure 
to that employed above, it is clear that one can determine functions of 
the form of £(— ki, . . . , — k v , k v+ %, . . . , k r \z) where kj > for all j < v, 
and kj > for v < j < r, but for more general tuples of both positive 
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and non-positive integers, the difference equation method of |Mes59] 
fails in that non-convergent series would arise as the supposed solu- 
tions. 



3.3. Algebraic independence of Hurwitz polyzeta functions. 

Theorem 3.11. The Hurwitz polyzeta functions are algebraically in- 
dependent over C - i.e for any N > 1, if P(Ti, . . . , T/v) G C[Xi, . . . , Tjy] 
is a polynomial which vanishes ( uniformly in z) at some N -tuple of any 
Hurwitz polyzeta functions of the variable z, then P is identically 
zero. 

Proof: Firstly, the linear independence of the Hurwitz polyzeta func- 
tions is an easy consequence of the difference equations Consider 
an arbitrary linear combination of a constant function and HURWITZ 
polyzeta functions of depths two and one, which sums to zero: 

A + ^2 A fclfc2 COi, k 2 ; z) + ^2 AfcCO; z) = 0. 

fci;fc2 k 

This equation is also valid replacing z by z + 1, so that from the func- 
tional equations ([3]), we find that also 

E A ^ 2 (-^f^ {z+ 1 m) k 2 ) -E A ^ = °- 

fei;fc 2 \ m=l v 1 / k 

But now it is clear that the coefficients A^^ and A^ must all be zero, 
since ^ has a pole of order exactly k at z = 0, and ^- Ylm=i ( z +m) fc 2 
has a pole of order &i at z = and a pole of order fc 2 at z = — 1. Then 
trivially, also A = 0. 

Now suppose that one has shown the linear independence of a C-linear 
combination of a constant function and HURWITZ polyzeta functions 
of depths 1, . . . , m — 1. For constants X,j p G C, consider now 



o = A+ V.JaC(iT,---,i 



m. \ 

mi " ) 



lm Jrt 

l-t .....L 



jm — 1 jm — 1 ^1 
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Because this relation must also hold for z+1, using (J3j) and subtracting 
the two equations from each other, we find that 

o= E ^•••c^C(^2V-- ; C; 2: + 1) 

]m Jm 

+ E \T~ l -i m ~\ im-iCtt™ S • • • > C-i; z + i) + 

r 1 m— 1 

jm— 1 jm-1 
'l '•••' t m-l 

Now let 

m 

'=nn<; 

and multiply the equation by z L . Then each function z L ~ l i((l 3 2j . . . , ij; z+ 

1) has a zero of order exactly L — l{ at z = 0. Functions with zeros 
of distinct orders at a point of C are linearly independent over C, so 
the equation breaks up into distinct parts where the ^-values agree. In 
this way a system of equations of the form: 

E v^'OT. • • • . c * + 1) + • • • + E M zL ~ l = 

/"i Im il 

results. (Here I denotes the common value of the indices I™, . . . 
But these equations can be divided by z L ~ l and then the inductive hy- 
pothesis guarantees that all of the coefficients are zero. This concludes 
the linear independence proof. 

Thanks to the stuffle product, the functions are also algebraically inde- 
pendent, since any supposed algebraic relation among such functions 
could be decomposed by means of the stuffle product into a linear ex- 
pression. In such a linear expression, the coefficients are sums of the 
coefficients of the original algebraic expression. However, at least one 
of the resulting polyzeta functions only arises from a single term of the 
original expression, the coefficient of which must therefore be zero. This 
eliminates certain terms from the linear expression, and once again, at 
least one of the remaining polyzeta functions comes out of a unique 
term of the algebraic expression, so that this latter term again has zero 
coefficient by the linear independence. Continuing inductively in this 
way, the theorem is proven. □ 



• • • + y^i-r. 
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